Image resonance in the many-body density of states at a metal surface 
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The electronic properties of a semi-infinite metal surface without a bulk gap are studied by a 
formalism able to account for the continuous spectrum of the system. The density of states at 
the surface is calculated within the GW approximation of many-body perturbation theory. We 
demonstrate the presence of an unoccupied surface resonance peaked at the position of the first 
image state. The resonance encompasses the whole Rydberg series of image states and cannot be 
resolved into individual peaks. Its origin is the shift in spectral weight when many-body correlation 
effects are taken into account. 

PACS numbers: 73.20.At 



I. INTRODUCTION 

For the understanding of several fundamental prop- 
erties of condensed matter surfaces, the knowledge of 
the electronic density of states and of the parallel- 
wavevector-resolved density are essential. The dispersion 
of both occupied and excited surface states, the influence 
of the changes of the electronic structure in thin film 
growth and enhanced magnetism are just a few exam- 
ples of relevant quantities. Experimentally the density 
of states can be probed by a variety of techniques such 
as photoemission spectroscopies^ (angle resolved, inte- 
grated, inverse and two-photon), and scanning tunneling 
microscopy^ To provide an adequate theoretical descrip- 
tion of the experimental observables, it is necessary to 
employ methods which retain the continuous character of 
the spectrum. In other words one has to take into account 
that certain quantities may be defined for any energy in a 
given energy interval. Furthermore, the formalism must 
be capable of describing excited-state properties. 

Regarding excited states, considerable experimental 
interest has been devoted to image-potential induced 
(IPI) states^*^ and resonances (I&iLiil IPI states are 
present in systems where a bulk band gap provides a 
barrier, trapping electrons in the image tail of the sur- 
face potential. If no gap is present at the IPI ener- 
gies, an electron is not reflected completely at the bulk 
barrier and hybridization with surface truncated bulk 
states becomes possible. This results in the forma- 
tion of resonances for some materials. A comprehensive 
theoretical description of IPI states has been given by 
Echenique and coworkers jiii^ii&ii whereas the situation 
for IPI resonances is much less satisfactory. In prin- 
ciple, IPI resonances are intrinsically contained in the 
many-body framework, already at the level of the GW 
approximation, 15 but to our knowledge surfaces have 
only been investigated in this context using a repeated 
slab geometryiA£iiLi£ Such a simplified treatment cannot 
capture the continuous spectrum of a real surface, be- 



cause the spectral function will inevitably be composed 
of a limited number of sharp, discrete peaks in place of 
the resonance. 

In this Article we present ab initio many-body calcu- 
lations of the local density of states (LDOS) of a semi- 
infinite jellium surface and demonstrate the presence of 
a broad IPI resonance. We calculate the LDOS - in its 
many-body generalization, the spectral function^ - de- 
composed according to the surface parallel wave-vector 
kii within the surface zy-plane as 



(1) 



where A is the k|| -resolved LDOS or spectral weight func- 
tion, defined by 



A{z,k»,u) 
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3G(z, z, k||, cj)sgn(o; — /i). (2) 



Here G is the one-particle Green's function in the rep- 
resentation indicated, and Hartree atomic units (ao — 
0.529 A, 1 hartree= 27.2 eV) are used. 
G is obtained from Dyson's equation 



G = G DFT + G DFT [S X c-^xc]G, 



(3) 



where «xc is the exchange and correlation DFT potential 
and Exc is the GW electron self-energy. 

Equation © is solved using a recent method developed 
to perform GW calculations in infinite, non-periodic 
geometries^ based on the embedding method, 2 » The ad- 
vantage of this approach is that the semi-infinite sub- 
strate, surface and vacuum regions are treated equally 
without the need for any fitting parameters or a repeated 
cell geometry. 

The main steps of the computation are as follows: 
(i) The Kohn-Sham equation^ is solved self-consistently 
within LDA. 22 In this framework, the embedding method 
permits an exact treatment of the semi-infinite substrate. 
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FIG. 1: Spectral weight function at the chemical potential 
and zero parallel wave-vector for semi-infinite jellium with 
r s — 2.07 ao. 

G DFT , obtained by numerical inversion of the Kohn- 
Sham Hamiltonian, describes a continuous spectrum of 
fictitious non-interacting electrons, (ii) From G DFT we 
compute the polarization P in the random phase approx- 
imation, (hi) The inverse dielectric response, (1 — vP)^ 1 , 
yields the effective interaction VF RPA . (iv) The self- 
energy Exc = iG DFT VF RPA is calculated using a real- 
space, imaginary-frequency representation and one ob- 
tains the self-energy on the real frequency axis by means 
of analytic continuation, (v) Equation @ can now be 
solved to update G. 

In Sec.[n]we present and discuss the results. Section lnTl 
is devoted to conclusions. 



II. RESULTS AND DISCUSSION 

We examine first the z-dependence of the spectral 
weight function A(z, k|| , u>) . This quantity is proportional 
to the probability amplitude for a particular wave- vector 
k|| and energy oj. In Fig. ^ we report A(z,\s.\\,uj) for 
ku equal to zero and u> equal to the chemical potential 
fi. All results shown in this Article are obtained for a 
jellium substrate of aluminum density (r s = 2.07 ao). 
The spectral weight outside the surface is enhanced by 
the improved description of exchange-correlation effects 
in GW, in common with the states of Al(lll) studied 
in Ref. l2.'il By varying the energy w we have verified 
that this feature is common to all bound states. In the 
bulk, the GW spectral weight is lower than that calcu- 
lated by the LDA because some weight is transferred to 
lower energies through electron-plasmon coupling. The 
larger amplitude of A in the surface layer, together with 
the absence of decay into the bulk, identifies the states 
in this part of the surface band structure as forming a 
surface resonance. 

In Fig. H we plot Act = a GW - ct lda , the difference 
between our many-body LDOS and that in the LDA. 
The energy-dependence of Act as z moves from the bulk 
towards the vacuum is indicated by contour levels. We 
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FIG. 2: Contour levels of Aa(z,u}) for semi-infinite jellium 
with r s = 2.07 ao. The solid horizontal line indicates z = 
3.5 ao (see Fig. 0J. The pronounced "teardrop" feature at 
z ~ 3.5 ao is the surface image resonance. Vv is the vacuum 
energy level, and Vb the bottom of the free-electron band. 



recapitulate that, for a jellium surface, the LDA effec- 
tive potential v e s(z) approaches the constant limits Vb 
and Vy for z — > — oo (bulk) and z — > +oo (vacuum) re- 
spectively. As a consequence, the bulk and vacuum lim- 
its of ct lda are simply proportional to the well-known 
expressions y/oj — Vb and yfuj — Vy- In the bulk, the 
LDA therefore predicts no states below Vb ■ The electron- 
plasmon coupling, automatically included in GW, is re- 
sponsible for moving states down in energy from above 
to below Vb, yielding the sub-band shown in Fig. for 
energies around 1 hartree below the chemical potential, 
as already demonstrated for the homogeneous electron 
gas^ 

The presence of the surface is still noticeable even some 
atomic units into the substrate through Friedel oscilla- 
tions, which are visible for bound energies from the con- 
tour levels of Fig. |21 and in the LDOS at a distance from 
the surface z = —10 ao in Fig. 

When z approaches the surface, coupling with plas- 
mons considerably reduces, eventually becoming negligi- 
ble a few atomic units outside the surface. As z moves 
into the vacuum, the electron density decays to zero. 
Self-energy effects are no longer present and the LDA 
description becomes exact. We then obtain Act = 0. 

Apart from surface-truncated bulk structures, Fig. [21 
also highlights a completely new feature in the form of a 
teardrop-shaped enhancement of the LDOS, localized in 
the near-surface region at energies between the chemical 
potential and the vacuum level. To highlight this feature 
we plot the energy dependence of Act for z = 3.5 ao 
in Fig. 21 (note that this corresponds to a cross section 
plot along a horizontal line through Fig. [2}. The peak 
in Act, located at about 1/32 hartree below the vacuum 
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FIG. 3: LDOS a GW and a LUA for semi-infinite jellium with 
r s — 2.07 ao at z = — 10 ao from the jellium edge. The 
"ripples" between Vv and Vb indicate the presence of Friedel 
oscillations. 



FIG. 5: Local density of states for a semi-infinite jellium with 
r s — 2.07 ao at z = 3.5 ao- The quantity reported in Fig.0]is 
the difference between the continuous and the dashed line. 
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FIG. 4: Difference in the LDOS Act for semi-infinite jellium 
with r„ = 2.07 ao at z = 3.5 ao from the surface. The dashed 
vertical line indicates the position of the first Rydberg state, 
1/32 hartree below the vacuum level. 



level (the position of the first Rydberg stateii), clearly 
identifies an IPI surface resonance. 

The presence of IPI resonances, even in the limiting 
case of a substrate without bulk reflectivity (the bulk 
reflectivity is associated with scattering by the atomic- 
cores 24 ), has been unclear. Simple model potentials in 
independent-particle approximations produced IPI res- 
onances for jellium substrates^ in some cases, but the 
presence of a clear resonance depended sensitively on the 
precise details of the chosen model. We remind that IPI 
states and resonances are a many-body phenomenon and 
can thus be included only in an ad hoc way in single- 
particle dynamic approximations. In contrast to these 
earlier, parameter-dependent results, our method is fully 
ab initio and includes many-body correlations. 

In the case of semi-infinite jellium, the zero bulk re- 
flectivity yields a large linewidth of the resonance, since 
hybridization with bulk states is not prevented at all. 
We estimate of the full width at half maximum (FWHM) 



from Act as 0.07 hartree (about 2 eV)^& We remark that 
in systems without a bulk band gap the hybridization 
with the continuous bulk band is the major contributor 
to the linewidth, as the self-energy alone accounts only 
for about one tenth of the total FWHMi4 

Owing to its large linewidth, the resonance shown in 
Fig. 0] encompasses the whole series of Rydberg image 
states^ which in the infinite bulk barrier model extends 
from Vy — 1/32 hartree to V\. The resonance is clearly 
visible in the difference plot between the many-body 
LDOS and the LDOS in the LDA. In the many-body 
LDOS itself, on the other hand, it resolves as a shoul- 
der feature just below the vacuum energy (Fig. [SJ rather 
than a distinct peak, because the LDA LDOS is already a 
strongly increasing function of energy in this range. The 
quantity reported in Fig.0|is the difference between two 
theoretical quantities and is therefore not directly acces- 
sible to experiment. However, it allows us to isolate an 
important contribution to the phenomenon - the many- 
body interaction - which could improve the description of 
the electronic structure of real semi-infinite surfaces be- 
yond the usual single-particle approaches such as Ref.l2rj 

We will now examine the origin of the resonance emerg- 
ing in Fig. 2] For this reason it is important to look at 
the spectral weight function A, whose integral over k|| 
is the LDOS [Eq. JJ)]. In Fig. we present the energy 
dependence of A(z, k|| , uS) for k|| = at the same position 
outside the jellium edge as Figs. 0] and |U(z = 3.5 ao). 

We first consider the LDA result. For non- interacting 
particles in a constant potential V the z-resolved spec- 
tral weight function is essentially one-dimensional and 
hence proportional to (u> — fc 2 /2 — V)~ x / 2 . Thus in the 
vacuum a singularity is present at Vv for k\\ = 0. The 
spectral weight then diminishes with increasing energy. 
But near the surface weight is transferred to bulk states 
decaying out into the vacuum: the spectral function has a 
peak very close to v e g(z), but singularities are no longer 
present. We emphasize that in a slab geometry Fig. 
would appear as a collection of delta functions, each delta 
corresponding to one of the discrete eigenstates. 
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FIG. 6: Spectral weight function for a semi-infinite jellium 
with r a — 2.07 ao, for ku = and z = 3.5 do- The integral 
over k|| of this function gives the LDOS reported in Fig- El 

When we describe the system in the interacting pic- 
ture (GW), we still observe bulk states that spill out 
into the vacuum. But we also have a new class of states, 
constituting the IPI resonance. They have substantial 
weight at and outside the surface. As a consequence 
spectral weight from the LDA states is transferred into 
these quasiparticle states. Spectral weight hence moves 
down from higher energies to the energies of the IPI res- 
onance. This produces the shoulder in A at Vy, and the 
displacement of the peak energy to lower energies. 

If fcii is increased from its zero value, the profile of A 
shown in Fig.[5]is shifted towards higher energies by /2, 
without being distorted too much (i.e., the dispersion of 
the states is nearly parabolic with effective mass equal 



to 1). When we integrate over ku to get the LDOS, one 
notices that for energies roughly below Vy the GW spec- 
tral weight A GW is greater than ^4 LDA for any value of 
fcii , leading to the positive Act shown in Fig.^J For higher 
energies, A GW — A LDA changes from negative to positive 
as fc|| increases. Negative contributions dominate and the 
corresponding Act is negative. 



III. CONCLUSIONS 

In summary, we have evaluated from first principles the 
local density of states of a semi-infinite simple metal sur- 
face. Many-body correlations are included via the GW 
self-energy. The substrate is described as truly semi- 
infinite, thus enabling a calculation of the continuous 
spectrum necessary to properly account for the hybridiza- 
tion between surface electronic states and bulk states. 
We demonstrate the presence of an IPI resonance just 
below the vacuum energy, encompassing the full series 
of image-states that would be present in a system with a 
surface bandgap. The origin of the resonance is explained 
in terms of the spectral weight tranfer to lower energies 
due to the inclusion of electron-electron correlation. 
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